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ABSTRACT 

Hysteresis is observed at second order phase transitions. Universal scaling for- 
mulae for the areas of hysteresis loops are written down. Critical exponents are 
defined, and related to other exponents for static and dynamic critical phenomena. 
These relations are verified with Langevin dynamics in both the critical and tricriti- 
cal mean-field models. A finite-size scaling relation is tested in the two-dimensional 
Ising model with heat-bath dynamics. 
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Hysteresis has been studied quantitatively for one hundred years now [1] . This 
phenomenon is well-known at first-order phase transitions. We find that it can also 
be observed at a second-order transition. Such critical hysteresis is characterised by 
certain scaling laws, i.e., universal exponents and scaling functions. The exponents 
are related to well-known static and dynamical critical exponents. Critical hysteresis 
is thus seen to be entirely consistent with the phenomenon of critical slowing down. 
Furthermore, when the critical point is a limit of a line of first order phase transitions, 
the scaling function for the hysteresis loop area is a smooth limit of those obtained 
along the line of transitions. 

The characteristic observation in hysteresis is that when a system is subjected 
to a cyclic external field, /i, its response, m, lags behind the field. Thiis the curve 
m{h) forms a loop. Although many different kinds of systems show hysteretic be- 
haviour, in this letter we shall use a language appropriate to magnetic systems. Thus, 
the external field will be a periodic magnetic field, with period T = 27r/c<;, and an 
amplitude ho. The response of the system will be characterised by an induced time- 
dependent magnetisation. All the statements we shall make have their analogues for 
other kinds of phase transitions, provided the external field is linearly coupled to the 
order parameter, which forms the response. 

Of interest in most experiments is the area of the hysteresis loop, defined as 



where the integral is over one period of the applied field. This is the energy dissipated 
by the system. Also of interest are the shapes of the hysteresis loops. There have been 
many early studies of the dependence of ^ on a; and ho [2] . Recent theoretical work 
has focussed on the low-temperature phase of 0{N) models[2,3]. In these models 
power laws have been derived for the scaling of hysteresis loop areas with ho and oj. 

The essential physics of hysteresis is the lag of the response behind the force. 
This is due to finite autocorrelation times. Whenever the inverse frequency of the 
external force is comparable to an appropriate autocorrelation time of the system, 
one expects to see hysteresis. Due to the phenomenon of critical slowing down, one 
must therefore sec hysteresis also at a second-order phase transition. We shall first 
outline the dynamics that we use. Next, assuming the existence of this phenomenon, 
we shall present scaling laws for the loop area and define critical exponents. Finally 
we shall present some examples of critical hysteresis. The scaling laws are seen to be 
satisfied with appropriate exponents. Further details will be presented elsewhere [4]. 

In this letter we shall take the dynamics to be given by a Langevin equation 
without conservation laws. We shall consider real one-component fields 0j on the 
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sites, z, of a d-dimensional lattice, whose interactions are specified by a Hamiltonian 
H. In terms of a Gaussian white noise Vi{t), with width F, the Langevin equation 
can be written as 

Note that this dynamics has no conserved quantities. We shall restrict our attention 
to the zero mode or magnetisation, (f> = "^^(pi- The argument is that this Fourier 
mode has the largest autocorrelation time, and therefore the long-time behaviour of 
the system is due to this. It is known [5,6] that all correlations functions arising 
from this one-mode dynamics are exactly equal to those obtained from the Euclidean 
quantum mechanics specified by the Euclidean action 



(3) 



where W stands for T5H/ Scj), and the prime denotes a derivative with respect to 0. 
The 'Planck's constant', is equal to T/L'^. The infinite volume limit of the theory is 
therefore obtained as the 'classical limit'. This is investigated by the Euler-Lagrange 
equations arising from the Euclidean action above, and reduces to the Langevin 
equation with the noise term dropped. This method has previously been used [6] to 
obtain the instanton solution in the dynamics. 

For technical reasons, we shall assume that the amplitude of the applied magnetic 
field is small and linear response theory is applicable. A fiucutuation-dissipation 
theorem [7] relates the energy dissipation to the the absorptive part of the complex 
susceptibility, x"{^)i through the relation 

A{ujM = \hlujx"{^)- (4) 

Along with the definition 

X"(a;) = Xm H dte-^^'{4>{t)m), (5) 

where the angular brackets denote averaging over the equilibrium ensemble, we now 
have the technical machinery to compute quantities of interest. Note that the average 
over the equilibrium ensemble in the formula above is a reflection of the fact that we 
are working in linear response theory. 

As an example, consider the most trivial of models — a Gaussian integral. The 
Hamiltonian is 

H = \ojlct>\ (6) 
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By a variety of methods one can construct the autocorrelation function in Langevin 
dynamics 

Cit) = e-*/^ (r = IK')- (7) 



This yields hysteresis loop areas 



A{u;, ho) 



1 hluT^ 

2 1 + a;2T2 ■ 



(8) 



The claim made earlier, that hysteresis vanishes only when the autocorrelation time 
does, is explicitly demonstrated in this model. Note also the 1/uj fall at large to. This 
behaviour has been observed, for example, in 0{N) models in dimensions d > 1 [3]. 
It will be seen again in critical hysteresis. Such behaviour is obtained whenever a 
pole in the spectral density determines the behaviour of A at large frequencies. 

Now we write a scaling relation at second-order phase transitions for the hys- 
teresis loop area as a function of oj and Hq. The necessary presence of an external 
magnetic field implies that we can tune 6 = {T — Tc) /Tc to zero, and let Hq control 
the approach to criticality. As a result, exponents corresponding to 9 do not appear 
in the scaling relations. We write 



In the limit ct; — 0, this must yield the usual relation between the magnetic field and 
the magnetisation. This implies 



Time-dependent phenomena are characterised by the scaling of the field autocorre- 
lation time T with increasing correlation length ^. The relation between these two 
quantities, r ~ defines a dynamical critical exponent z. As one approaches the 
critical point by tuning /iq, the scaling of uj required to keep the hysteresis loop area 
fixed must be given by the relation between r and ^ above. This gives the remaining 
critical exponent 

The scaled variable ujHq shall be called VL. For non-zero 9, there will be, of course, 
corrections to the scaling form given above which become apparent as /iq ^ 0. 

A trivial illustration of the scaling formula is provided by the Gaussian integral 
analysed before. Coupling a magnetic field linearly to it is easy to see that 5=1 
and, since ujq is independent of /i, v/5(3 = 0. As a result, A = 2, consistent with the 



A{uj,ho,9) = 



(9) 




(10) 
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overall factor of h^, and = 0, consistent with the fact that the only ho dependence 
is in this factor. 

The scaling laws illuminate the origin of critical hysteresis. Provided that the 
scaling function F{0,) is non-zero, one must have critical hysteresis. This is a con- 
sequence of critical slowing down. We anticipate the subsequent discussion leading 
up to the fact that F{Q) has a maximum at (say) Qq, decays for Q > Qq as 
and for small Q increases as a power of Q (see Fig. 1). As a consequence, the scaling 
form above predicts that as ho decreases, the hysteresis loop area increases at any 
fixed arbitrarily small frequency. At large frequencies, on the other hand, the system 
cannot follow changes in the driving field, and the magnetisation remains zero. In 
this qualitative sense, hysteretic behaviour at a critical point closely resembles that 
at a first order transition. 

The scaling relations for a thermodynamically large system imply, in the usual 
way, finite-size scaling (FSS) relations. Well-known arguments [8] can be adapted to 
obtain various scaling formulae. Of particular interest is the finite-size scaling of the 
frequency, cu™,, at which the hysteresis loop area reaches its maximum. When systems 
of different sizes, L, are studied, each at some appropriately defined pseudo-critical 
coupling, and Hq is independent of L but smaller than the value at which finite-size 
rounding sets in for the largest system, then 



Note that this relation gives a method for the measurement of z. The usual FSS 
methods for dynamics only measure the combination of exponents z/u. 

Our first illustrative example is of the time-dependent Landau- Ginzburg model 
at its critical point. The model is defined by the Langevin equation in Eq. (2) along 
with the zero-mode Hamiltonian 



A (possibly time-dependent) magnetic field h has been introduced. The critical 
theory corresponds to r = 0. We investigate this theory for large volume systems. 
As explained earlier, this requires the solution of the ordinary differential equation 
obtained by dropping the noise term in the Langevin equation. After the scalings 



(jJm{L) ~ L^. 



(12) 



(13) 





(14) 



we obtain the equation 



dcf) 
Iti 



(j)^ -ce(j) + h{t). 



(15) 
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We used a harmonic form, ho cosut, for the driving field h{t). Note that the scahng of 
t imphes one for the frequency uj. We have assumed that r 0, and that c includes a 
factor y^GT/g, whose temperature-dependence can be neglected for sufficiently small 
9. The presence of hysteresis makes this a stiff equation. It was solved using Gear's 
method, starting from an arbitrary initial value and integrating until transients died 
out to the precision of the arithmetic. The subsequent time history of the system gave 
the stable limit cycles corresponding to hysteresis. Scaling was checked by varying 
ho over three decades and lv over four. 

The solutions showed hysteretic behaviour quite clearly. At very high frequen- 
cies, the system could not follow the applied magnetic field; and the magnetisation 
remained zero and time independent. As co decreased, the hysteresis loop m{h) 
opened into an ellipse with its major axis along h. With further decrease in uj, the 
major axis turned, and the ellipse began to get dented, until it resembled familiar 
hysteresis loops. Eventually, with decreasing uj these loops became thinner. We cal- 
culated the areas of the hysteresis loops as a function of uj and Hq and found that the 
scaling law given in Eq. (9) is obtained with high precision. We found the exponents 

4 2 
A=-, and (1^) 

These are completely consistent with the values S = S, u = jS = 1/2 and z = 2. The 
scaling function can be identified by our computation. It is shown in Fig. 1. It has 
a single maximum at Q ~ 0.93. 

These results were obtained for cd = 0. Positive values for this parameter 
correspond to approaching the critical point from the disordered side. As the value 
of this parameter increased, distortions in F{Q,) could be seen at higher and higher 
Cl. For small ft these tended to decrease F, and push flo towards larger values. 
For c9 less than zero, we approach criticality from the ordered side, where there are 
coexisting phases. Varying c9 in this region, for any fixed ho, we found that A{u!, ho) 
approached its value at criticality smoothly. 

It is interesting that the hysteresis loops themselves are scaling functions. At 
^ = it was found that 

m{h;u},ho) = hQ^^ L{h/ho',^), (17) 

where the loop shape L depended on Q. This is a form of a dynamical principle of 
corresponding states. The shapes of the loops are shown in Fig. 1 at the appropriate 
points below the scaled frequency axis. For — > oo, L = identically, and for — > 0, 
L{x) = x^l^. 
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The exponents and scaling functions enjoy a high degree of universality. We 
checked this by changing the form of the cyclic driving field h{t). For a variety of 
forms, all with time-reversal invariance, we found the same values for the exponents. 
Whether we used triangular waves, or piecewise constant functions did not matter. 
Furthermore, the scaling function obtained for each shape of the driving field was 
simply related to that obtained for the harmonic form. We found that all the scaling 
functions mapped on to each other with multiplicative redefinitions of ho and uj. 
This leads to the form ciF(c2fi), with universal F{0,), where ci and C2 are constants 
independent of ho, u> and 9, dependent only on the shape of the driving force. 

A second example is provided by the Langevin dynamics of the (f)^ model at its 
tricritical point. After appropriate scalings, we have, for the large-volume limit of 
the Langevin equation 

^ = -cf>' + h{t). (18) 

Numerical solutions were obtained for this equation, varying Hq and u as before. 
The critical exponents are well-known [9] — S — /3 = 1/4, v = 1/2 and z = 2. 
These yield A = 6/5 and /x = —4/5. These described the numerical results extremely 
well. The universal scaling function thus obtained is shown in Fig. 1. As for the 0^ 
model, this function is universal over cyclic driving fields. It is also interesting that 
the scaling functions for critical and tricritical dynamics are identical for Q ^ Q^- 
Universal hysteresis loop shapes were also observed, and found to be qualitatively 
similiar to those for the 0^ theory. 

We checked the FSS relation in Eq. (12) through heat-bath simulations of the 
two-dimensional Ising model. We simulated 32^, 40^ and 64? lattices at the pseudo- 
critical coupling, P^, defined by the maximum of the specific heat. Equlibration at 
zero field was achieved by 10^ sweeps with the Swendson-Wang cluster algorithm. 
Hysteresis runs were performed changing the external field cyclically between ±ho 
in steps of Ah. On all the lattices we chose ho = 8 x 10~^ and Ah = 10~^. At 
each value of the field we ran N steps of the heat-bath algorithm. The frequency 
is given by a; = 27rAh/Nho. Expectation values and errors for hysteresis loop areas 
were obtained by jack-knifing measurements from 500 loop traversals at each lv and 
L into 10 blocks. A series of runs on each size of lattice located tOm with fairly 
high precision. The data are shown in in Fig. 2 along with the scaling law obtained 
using the dynamical exponent z = 2 for heat-bath dynamics. The hysteresis loop 
shapes were seen to be universal, and qualitatively similiar to the time-dependent 
Landau- Ginzburg model. 

In conclusion, we have observed hysteresis at a critical point. Scaling relations 
and critical exponents for hysteresis loop areas have been identified and related to 



6 



other static and dynamic critical exponents. Finite-size scaling relations following 
from these have been written down. Since the observation of hysteresis is no guarantee 
that a first-order phase transition is seen, these relations are important in deciding the 
nature of the transition. Scaled hysteresis loop shapes and areas are also universal, 
characterised by functions which do not depend on the precise form of the driving 
force. The scaling functions identified here should be interesting objects of study, 
since they contain information on the complex structure of the underlying spectral 
density in an easily accessible fashion. 
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The universal scaling function = ^((^/iq) describing hysteresis in the 

time-dependent Landau-Ginzburg theory. The result for the critical 0^ theory is 
shown with a full line and that for the tricritical 0^ theory with a dashed line. The 
universal loop shapes L{h/ho',^) in the model are shown below the main figure. 
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FIGURE 2 

Finite-size scaling tests for hysteresis in the two-dimensional Ising model. The 
data have been obtained at couplings corresponding to the maximum of the specific 
heat. The line shows the FSS prediction using the dynamical exponent z = 2. 
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